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Abstract. 

All models of interacting electrons and spins can be reformulated as theories of 
interacting Majorana fermions. We consider the Kondo lattice model that admits a 
synirnctric representation in terms of Majorana fermions. In the first part of this work 
we study two variational states, which are natural in the Majorana formulation. At 
weak coupling a state in which three Majorana fermions tend to propagate together 
as bound objects is favored, while for strong coupling a better description is obtained 
by having deconfincd Majorana fermions. This way of looking at the Kondo lattice 
offers an alternative phenomenological description of this model. In the second part of 
the paper we provide a detailed derivation of the discretized path integral formulation 
of any Majorana fermion theory. This general formulation will be useful as a starting 
point for further studies, such as Quantum Monte Carlo, perturbative expansions, and 
Renormalization Group analysis. As an example we use this path integral formalism 
to formulate a finite temperature variational calculation, which generalizes the ground 
state variational calculation of the first part. This calculation shows how the formation 
of three-body bound states of Majorana fermions can be handled in the path integral 
formalism. 
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1. Introduction 

One of the most studied models in condensed matter physics is the Kondo lattice model, 
which consists of a lattice of localized magnetic moments (/-spins) interacting with a 
sea of conduction electrons (c-electrons), via an exchange coupling. This situation is 
characteristic of many varieties of rare earth and actinide compounds, and it is believed 
to be responsible for the incredibly rich physics of these systems. The properties of this 
model system have been studied intensively both analytically and numerically, deploying 
many different techniques, see e.g. the reviews [1, 2, 3, 4]. The main issue encountered 
by any analysis of the Kondo lattice model is the difference between the algebras of the 
/-spins and the c-electrons. Because of this difference it is very hard to describe the 
two subsystems on equal footing. 

In a recent work [5] we introduced a faithful fermionic representation of the 
Kondo lattice, that was based on the representation of the Hamiltonian in terms 
of Majorana fermion degrees of freedom (Majoranas). Using this representation it 
is possible to treat spins and electrons in a more symmetric way, without focusing 
on one subsystem in particular. The common (Clifford) algebra fulfilled by the 
Majorana fermions allows for the generation of electron-spin canonical transformations, 
which are otherwise quite involved when written in terms of the standard operators. 
Majorana fermion representations of spins is an old technique [6, 7] that was recently 
popularized by its application in the influential Kitaev model [8]. In the context of our 
study this representation was used in works focusing on non-Fermi liquid behavior in 
modified Kondo impurity problems and lattice systems [9, 10, 11]. The fact that this 
representation can be used to study spin models and generate spin liquid states has also 
been known for quite some time [12, 13, 14, 15]. To the best of our knowledge [5] is the 
first application of Majoranas in the standard Kondo lattice model. 

In this article we are going to improve on our previous results in two ways. In a 
first stage we will introduce a canonical transformation among the Majorana fermions, 
that will allow us to obtain better results for the upper bound on the ground state 
energy. To do this we consider the limits of weak and strong coupling. We will work 
under the assumption that at weak coupling the three Majorana fermions coming from 
the /-spins, stick together to form a coherent Majorana fermion, that becomes (to a 
good approximation) the appropriate degree of freedom relevant in this limit. This 
interpretation is naturally suggested by the form of the Hamiltonian. Increasing the 
coupling strength, the coherence among the three Majoranas decreases, and hence it 
becomes more natural to think of them as independent particles. We indicate this 
process with the term deconfinement of the Majorana fermions. Both these descriptions 
are exact in the extreme limits of zero coupling and no hopping (atomic limit). In the 
second part of the paper we provide a detailed derivation of the imaginary time path 
integral formulation for a general Majorana fermion Hamiltonian. This will permit us 
to study the Kondo lattice model at non-zero temperatures and to take into account 
entropic effects. To do this we will discuss in detail how to handle the three-body 
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Majorana bound states, taking into account the effects due to the deconfinement process 
at a mean field level. 

The paper is organized as follows: in section 2 we briefly review the Majorana 
representation of the Kondo lattice; in section 3 we will then immediately introduce 
the canonical transformation among the Majorana fermions and use it to find the best 
variational guess for the ground state. In section 4 we study the path integral formulation 
of a general theory expressed in terms of Majorana fermions and in section 5 we will 
make use of this formalism to study our Kondo lattice Hamiltonian. At the end, in 
section 6, we provide a brief outlook and a summary of the results. 

2. Majorana fermion representation of the Kondo lattice model 

The Hamiltonian for the Kondo lattice, that we are going to discuss in this paper, is 
conventionally written as [1] 

Hklm = - ^ XI 5Z 4,<7(ri)cc,a(ri) + h-c (1) 
- [nc(rj) - 1] + J X Sc(ri) ■ S/(ri). 

i i 

The first line describes the hopping of the conduction electrons on the lattice, and a =t, J, 
is the spin label. In the second line the last term represents the (on-site) interaction 
between the local magnetic moment and the spin of the conduction electron, while the 
first term is the chemical potential term for the conduction electrons, so that nc(rj) = 
Xlo- o-(''i)cc,o-(rj) is the number operator. In the following we will consider only the half- 
filled Kondo lattice model, that implies a value of zero for the chemical potential. The 
operator S/(rj) represents the spin-| object coming from the magnetic moment on site r^, 
that satisfies the usual SU(2) algebra [S''{vi), S\rj)\ = idije^^^S'^ivi) and = 3/4 [16]. 
Here, and in the following, we will work in units such that h = 1. The c-electron's spin 
operator is given by the standard representation Sc(rj) = ^ Xlo-o-' 4,(T(i"j)'''o-,o-'Cc,cr'(rj), 
with r° (a = 1,2,3) the Pauli matrices [16]. 

Following [5] it is possible to represent the Hamiltonian (1) in terms of Majorana 
fermions. In order to do that, it is necessary to decompose the c-electron and the 
/-spins in terms of Majorana fermion degrees of freedom. It is well known (see, e.g. 
[10, 13]), that to represent spin-^ operators, three species of Majorana fermions are 
needed: /ia(rj), (a = 1,2,3). The Majorana fermions are real IJ,\ = Ha and independent 
{/^a(i'i)5 /^b('^i)} = ^ij^ab- It's Straightforward to check that the operators 

5-(r.) = -^e<^^>,(r,)/x,(r,)/2, (2) 

satisfy the angular momentum algebra for spin-| operators. It is fundamental to note 
that a fourth Majorana fermion can be formed as a non-linear combination of the three 
just introduced: 

7o(ri) = 2i/ii(ri)/i2(ri)/i3(ri). (3) 
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This is a proper Majorana fermion (with 7q = 1/2), that commutes with the /-spin 
operator [7o(rj), S/(rj)] = 0. As in [5] we can now use this Majorana fermion, together 
with three other Majorana fermions, to represent the c-fermion operators: 

c,(r.) = e-^IiM^pM. c,(r.)=e-? -^''---^'^°'---> . (4) 

The difference from the usual Majorana representation of conventional fermion operators 
[10, 17], lies in the fact that 7o(rj) is here the composite operator (3). The exponential 
factors in (4) generate a site-dependent phase. This phase convention can be generalized 
to other bipartite lattices, it is here written for a simple cubic lattice and is useful to 
simplify the algebra and to emphasize the interpretation of 70 as the appropriate degree 
of freedom at small values of the coupling constant. After straightforward algebra, 
making use of (2)-(4) and assuming nearest-neighbor hopping on a simple cubic lattice, 
the Hamiltonian (1) can be rewritten as 



Hklm = H^'^ + H^'^ + H^'\ (5) 



with 



3D ^3 
a=l i d=l a=l i 

H^'^^ = - ^ ^[i7i(ri)/ii(ri)i72(ri);U2(ri) + permutations 1 ^ 2 ^ 3], (7) 

i 

where represents the lattice vectors (one for each dimension). A part in the kinetic 
term of the c-electrons is quadratic in 70 's, and it is therefore an operator of sixth order 
in iia- 

D 

^^'^=^^EE^o(r. + x,)7o(r.) (8) 

i d=l 

D 

= - z4t^^/ii(rj + Xrf)/i2(ri + Xrf)/i3(ri + Xd)/ii(ri)/X2(ri)/^3(ri). 

i d=l 

This formula is a direct consequence of the gauge choice (4) and points out the nature 
of 7o as a composite propagating object: the coherent nearest-neighbor hopping is the 
effect that dominates the physics of the /Za-Majoranas at weak coupling J/t. Therefore 
in this limit it is adequate to consider 70 and not the three independent /la-Majoranas 
as the appropriate degree of freedom. Vice- versa, in the strong couphng limit, the fia- 
Majoranas shall be considered as independent objects that hybridize with the 7a's, and 
the operator in (8) treated as a weak gluing interaction. It is important to note that 
both these states do not break the spin rotational symmetry, that is present in the initial 
Kondo lattice Hamiltonian, while in [5] the only good trial wave functions in the small 
J/t limit that we considered where broken symmetry states. To improve our results, 
keeping spin rotational symmetry, we will consider the hybridization of 7 and /i degrees 
of freedom in both limits. To optimize the linear combination between the two species 
of Majoranas, we will make use of a variational mean-field approach in the rest of the 
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paper, though the path integral formahsm introduced in section 4 permits to go beyond 
it. 

For future convenience, we now introduce the conventions that we will use to 
describe the system in momentum space. A generic Majorana fermion can be expressed 
in the momentum b 

7«(rO = ;^E'[^^''"7a(k) + e-^'^'-TKk)] = ^ E^^''"7a(k). (9) 

The prime in the sum denotes that only one of k and — k should be included in the 
sum. This convention brings to more readable formulas, because the operators 7(k) 
and 7^(k) fulfill the algebra of standard fermionic operators {7a(k), 7^(k')} = 5kk"^a6 
and {7a(k), 7f,(k')} = 0. However, it is often convenient to use the last expression 
on the r.h.s. in the intermediate steps. We will, for simplicity, consider lattices 
with Nd = 2Mrf and odd in the following. Then the allowed values of k are 
k = 27r( "^^^^^ , "^TVa'^^ '" • •)• This means that k and — k are distinct and that 

zero-modes are absent. The operators that appear in (6) can readily be reformulated in 
k-space, via a Fourier transform. The first term reads 

3D 3 ^ / 

^^EEE7'^(^^ + ^'^)7'^(^^) = i^o + EE'^k7l(k)7a(k), (10) 

a=l i d=l a=l k 

were Eq is the ground state energy for J = [in ID Eq/N = —At/n in the 
thermodynamic limit oo) and 

ek = 2t^sin(xd-k) = 2t J]sin(/crf). (11) 

d d 

It is natural to choose the k's in the sum so that ek > and hence the ground state has no 
7-excitations. For example in ID this implies that < < vr, while in 2D it corresponds 
to the square enclosed by the lines —k^ < ky < 27i — and — vr + k^ < ky < n + k^- 
We will identify these subsets of the full Brillouin zone with the symbol BZ'. The other 
quadratic term in the Hamiltonian is the hybridization term, which assumes a simple 
form in momentum space: 

^EE^7a(r.)/i.(r.) = ^X^$^'[.7l(k)/i.(k) -^/it(k)^^(k)]. (12) 

a=l i a=l k 

The operators that appear in (7) and (8) can also be expressed in k-space. However, 
in this case it is more convenient to postpone their treatment until we have defined 
our mean-field decomposition scheme, to simplify these interaction terms as much as 
possible. 



3. Variational calculation 



To gain energy from the quadratic part of the Hamiltonian H^"^^ [in particular from (12)] 
in the Majorana basis, it is natural to perform a rotation between 7's and /i's. To do 
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this we define rotated operators (a = 1, 2, 3) for eacli value of k in BZ': 

/ia(k) = cos(ak/2)/i„(k) + i sin(ak/2)7^(k), ^^^^ 

7a(k) = cos(ak/2)7a(k) +«sin(ak/2)/i„(k). 

Note tliat tfiis implies that a;_k = — Ok since we want to preserve standard anti- 
commutation relations. Applying the transformation to the quadratic piece we obtain 

H^"^ = EE'[^~7(k)7l(k)7,(k) + 6~,(k)/it(k)/i,(k)] 

a=l k 

3 

+ E E'^W[z7t(k)/i„(k) - ^/il(k)7„(k)] + 3Eo/4, (14) 

a=l k 

where for generic ak we have 

e^(k) = cos^(ak/2)ek + sin(ak)<//4, 

e^(k) = sin2(ak/2)ek - sin(ak) J/4, (15) 
y(k) = cos(ak)<>'^/4 — sin(ak)ek/2. 

Note that e/i(k) < when < tan(Q;k/2) < J/(2ek). The choice tan(ak) = -^/(2ek) 
diagonalizes the quadratic part of the Hamiltonian, i.e., it sets V^(k) = 0. However, as 
we shall see, this is not a good choice for small J/t since the cost in kinetic energy of the 
7o's is then too large. It is also convenient to introduce real space Majorana operators 
corresponding to the /i(k)'s and the 7(k)'s. To do this we use the same convention as in 
(9), but with tildes on the operators. We also define rotated versions of the three-body 
bound state and spins via 

7o(ri) = 22/ii(ri)/i2(ri)/i3(ri), 5'„(ri) = -ieabcfJ'b{ri)fic{ri)/2. (16) 

In the real space basis the relation between operators are (a = 1, 2, 3) 



(17) 



3 

with matrix elements given by 

2 / 2 / 

Aij = ^E cos(ak/2) cos(k-rjj), Bij = sin(ak/2) sin(k-rij).(18) 

k k 

Since the Hamiltonian is invariant under inversion k — t- tt — k it is natural to assume 
that Ok = tt7r-k, which implies unbroken inversion symmetry. A consequence of this 
choice is that the matrix elements of A (B) are nonzero only when the two indexes 
belong to the same (different) sublattice. This observation will simplify the algebra in 
the following. 
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3.1. Trial state with three-body bound states and totally uncorrelated spins 

We will consider two simple trial states for the ground state of the half-filled ID Kondo 
lattice, although similar calculations are easily performed in higher dimensions as well. 
Both states do not break spin rotational symmetry. The first trial state, which works 
better for small J/t, has a definite fermion parity for each of the flavors 7^ (a = 0, 1, 2, 3). 
This implies that all trial state averages involving an odd number of 7a's will vanish. The 
state we consider has a definite occupation of each fermion state and totally uncorrelated 
spins, explicitly 

(n^,(k)) =Oorl, for a = 0,1, 2, 3, 

(5„(r,)^,(r,)) = «,/4. 

The variational energy of this state is easily expressed in terms of the A and B matrices 
in real space and the angles ctk- The averages we need are 

3 

{H^"^) = E E'G~7(k)(n,.(k)) + 6,(k)/2) + 3E0/4, (20) 

0=1 k 

{H^'^) = - ^(E'«i^(«OK(k) - I)) (E'sin(ak)K(k) - i)) 

k k 

+ cyclic permutations 1 ^ 2 ^ 3 ^ 1, (21) 
(^(^)) = t J2 AfjAlY! 2sm{k ■ (r,, + f,))K(k) - 1/2) 

i,j,l,d k 

3 

- 4t7VEn(]^E2«in(Msin2K/2)(n,„(k) - 1/2)). (22) 

d a=l k 

We now consider the state with zero fermion excitations {n^^{\i)) — 0, since it has the 
lowest energy. The resulting energy functional (measured with respect to the ground 
state at J = 0) becomes 

k 

3J / 1 v-^' . , ,\2 At ,0 
- y(]7E + ^^00^02 + ■■■ (23) 

k 

In principle the full expression should be optimized with respect to ctk, but a detailed 
analysis shows that for small values of J/t the result is dominated by the first two terms. 
Moreover, in this limit we find that Aij,Bij,<^ Aqo ~ 1 for all i ^ j . Minimizing 
just those two terms, using the fact that ^00 ~ 1, the optimal angles are to a good 
approximation given by the solution to 

— sin(Q!k/2) + ^ sin(Q;k) = y cos(Q!k). (24) 

An even simpler trial state is to take a constant angle ctk = a. In this case the total 
energy functional is 

sin^fa) 

= 3 — sm [a/2) — — sm(Q;) — h 



N Vtt ' ' ' 16 ' ' 32 
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+ -[1 - cos^(a/2)] + % sin^(a/2), (25) 

and the optimization can be performed analytically. In figure 1 we plot the variational 
energy with numerically optimized c^k for the full Ai?t,ar[ak]; which is slightly better than 
the restricted variational states corresponding to (24) and (25). It is also compared with 
the conventional antiferromagnetic state in which the /-spins are locked into a static Neel 
order. In contrast to the trial states that we consider this state breaks spin rotational 
invariance. The variational energy of this state is given by [18] 

^ = -^E'(\/(W + ^W-^(k)). (26) 

k 

Both the energy of the Neel state and the the energy of our trial state (23) reproduce 
the correct leading constant term in the limit J/t — )■ 0. The subleading corrections make 
the Neel state slightly lower in energy for small but finite J /t. 

3.2. Trial state with independent deconfined Majorana excitations 

For large values of J/t it is important to gain energy from the dominating local term. 
This can be achieved by considering the trial state with 

(n^„(k))=0 for a = 1,2, 3, 

(n^^(k)) = l for a = 1,2, 3. 

The calculation of the trial energy functional of this state is straightforward, with the 
result 

k k 

- 4* E (^E' ^i^^^'^) cos(ak)) (28) 

d k 

This is the 0(3) state introduced in [5] in a different notation. It is also displayed in 
figure 1 and is clearly favored with respect to the other considered states for large values 
of J/t. In fact it gives the leading term correctly in the limit J/t ^ oo. 

3.3. Trial states with non-trivial spin correlations 

The trial state of section 3.1 can be further improved by allowing for non-trivial spin 
correlations. Indeed, if we consider a state where the fermion subsystem is in the 
ground state, but leaving the state of the (rotated) spin system undetermined, the 
expectation value of the Hamiltonian will generate a spin Hamiltonian that depends 
on the variational parameters ak- The best trial state is then obtained by setting 
the spin state to the ground state of this spin Hamiltonian. Optimizing the resulting 
trial energy with respect to the best possible trial state in this class is obtained. 
We leave this calculation for a later study since finding the ground state of the spin 
Hamiltonian is a highly non-trivial problem. We have checked that antiferromagnetic 
spin-spin correlations will improve the variational energy though. 
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Figure 1. Variational energies for the half-filled ID Kondo lattice for the different 
variational states. Solid red line: deconfined Majorana state; Dashed blue line: state 
with confined 3-body Majorana bound state; Dotted yellow line: broken symmetry 
Neel ordered state. 



4. Path integral formulation for Majorana fermions 

Using standard path integral formulas and fermion coherent states [19] we can calculate 
the expectation value for an operator O given the density matrix p = e~^^ /Zq from 

X {iMWi){-Ud\i,). (29) 

Zq = e'^° is a constant that is needed to insure that Trp = 1. The formalism can also 
easily be adapted to the calculation of the partition function, in which case we substitute 
p e"^^ and set = 1. Now we write p = p^/^ p^/^ ■ ■ ■ p^/^ and insert M - 1 identity 
operators with the slightly non-standard convention 

1 = / (n^C/C,m)e"^"^^''^'^^'^H^w)(6|. (30) 
The standard (Trotter) procedure then gives (Ar = (3/M) 

= Q-T,cil,Mio,,ie~^^^*'^\ (31) 

up to corrections of order 1/M. Here we have also used the discrete Fourier convention 
(with / = 1,...,M) 

The j-sum is over the M distinct solutions to e*'^-'*^ = —1, corresponding to anti- 
periodic boundary conditions. We will take M to be even so that all positive 
frequencies can be paired up with a corresponding negative one U-j = —ojj. This is 
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a unitary transformation and hence the corresponding Jacobian is just a phase. The 
proportionahty sign in (31) is needed since we have left out the integration over the 
internal variables. After rearrangement the full integration measure becomes 



n <,o^C,o) (n di.,.,dC,.) , (33) 



and the path integral representation for the trace is 

Ttdp = J d[^]e-^-(«-.A/+C,o)(€.,i+«a,o)c)(_^*^^je-^(«*'«)e-^o. (34) 

A matrix H that allows for a frequency dependent self-energy is obtained via the 
substitution 

E ^-^Hi^l^,) |. E E (35) 

I a j 

whereas for a diagonal quadratic Hamiltonian we will have 60,(0;^) = e^, independent of 
Uj. Introducing the propagator via the relation G~^{uj) = —[1 — e"""^^' — f3ea{oL!j)/M] 
we may then write the quadratic part of the action as 

5(r,0 = EC^.^«'(^^)^",-r (36) 

This implies that the action is diagonal in the a-index, i.e., we can write p = Yla Pa- AH 
of these are well-known results, let us now adapt this formalism to Majorana operators. 

4.I. Majorana path integral 

Let us now consider a Hamiltonian that is defined in terms of a set of Majorana 
fermion operators Xa- To use standard fermion coherent states the Majoranas have to 
be paired up to form conventional Dirac fermions. There are two straightforward ways 
to do this, called "fermion doubling" and "fermion halving" in [20]. In fermion halving 
the Xa themselves are paired up, this procedure is outlined in [21]. We use fermion 
doubling, which is more symmetric and relies on an extension of the Hilbert space. For 
each Xa we introduce a conventional fermion annihilation operator and represent Xa 
by :^(c^ + cjj). The other Majorana operators -^{c^ — c]^ do not enter the Hamiltonian 
and can be paired up to form an independent Hilbert space that will generate an overall 
degeneracy factor of 2^^/^ in the trace. To avoid fractional degeneracy we will take A^^ 
to be even, this is also always the case in physical applications where Majorana fermions 
unconditionally appear in pairs. This construction means that we have artificially 
doubled the size of the Hilbert space; we will eliminate these extra degrees of freedom 
from the path integral later. We can now directly use the path integral expressions 
above and perform a basis change for the Grassmann numbers to 

^ i + ^* / i i~ C I 

^ — \/2 ' ' ^ — \/2 ' ' d^a,id^a,i = di'a,idCa,i, (37) 
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The integrand of (34) then goes into 

Expanding the new Grassmann variables as Ca,i = -ttj Tlij ^~^^^^Ca,u}j (and the same for 



Va,i) the action becomes 

•5(C, ^) = Y. + - cos(a;,-)](^a,--,Ca,-, - Ca, 

+ ^^isin(u;j)(i/„_^^z/a 

a J 

Here the prime on the sum means that only positive frequencies are to be included. 
Expanding out a generic 0(Co) we have (repeated indexes are summed over) 

C(Co) = CafiH^al'^a' fl + H^ap-ysCaflCpfiC^iflCsfl 

+ H^'^l^^'pi'^'CaflCpflC'yflCa'flCp'fiC'y'fi + • • • (41) 

The integral over variables and ^* q for a's that are not present in the expansion of 
(9(Co) gives one together with the prefactor, i.e., 

t^Co^^Coe"^^"''^^^^"'"^^^"-^^^"-"^ = 1- (42) 
For a's that are present in 0(Co) the prefactor can be replaced by 

Expanding this and performing the integral over (^0,0 and Vafl imphes that we can 
substitute 



Ca,0 



V2 

1 ^ ^-^u,/2 [cos(a;,/2)C„,., - i sin(a;,/2)z/,,^J , (44) 



M 

3 



in 0(Co)- This implies that the integrand goes into 

O ( ^"■'""^"■M e-^K'^^^e-^" . (45) 



The corresponding integral (average) is easily evaluated by taking derivatives of the 
generating function with respect to the Grassmann source fields 77^: 

rf[e]e-^(^'")e^'^^" (46) 

where the measure can be written as 

d[i\ = n n '^^<^-^j'^^c,,uj,dCa-ui^dCa,oj^ = d[iy]d[(]. (47) 



a J 
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Because the exponent is now a quadratic polynomial in the z/'s, these variables can 
be integrated out exactly. Therefore the generating function is given in terms of an 
integral over Grassmann variables obtained by replacing, locally in the action, the 
Majorana operators in the Hamiltonian with Grassmann numbers. The imaginary time 
evolution of these numbers deviate from those of conventional fermions. In fact, after 
some straightforward algebra, the result is (dropping the overall normalization factor) 

Z(,,) = j dlCle-^K-'), (48) 

s(c. = A. 5: i/(c,) - . (y) C"^-"/-, + E %T.r^.- 

I j,a a ^ j 

This result can be simplified further by a linear transformation to another set of 
Grassmann numbers defined via 

-iWj/2A 

y = - ^ f49) 

^ COS[UJj/2) 

or equivalently, in discrete imaginary time 

Ca,l = ^ [Xa,l + (1 - Sl^M)Xa,l+l " Sl^MXa,l] ■ (50) 

In terms of these variables the generating function is 

Z{v) = f d[x]e-'^''''^\ (51) 



2 —-"^^-""-^ ^ i]^ 

In this expression it is clear that the source term is localized exactly at time zero. 
The main result of this section is that it is possible to reduce the calculation of the 
partition function (or a density matrix) for a Hamiltonian that is expressed in terms 
of Majorana fermions, to a standard Grassmann integral obtained by substituting 
the Majorana operators Xa with Grassmann numbers Xa,i for each time slice /. The 
properly regularized (Trotter) generating function is then given by (48) or (51). The 
corresponding continuum short-hand expression is (cf. [14, 21]) 

S{X,V)= f dr(^^J2x-^''^^rX-{r) + H[x{T)] + J2va{T)xa{r)), (52) 

a a 

where we have also included source fields at intermediate times. The regularization in 
(51) is a kind of "midpoint rule" [20], which is clear from (50). For the calculation of the 
actual value of the partition function a prefactor equal to 2^^^^^'^^^"^ should be included 
on the right hand side of (51). This can be derived by considering the partition function 
for H = T]a = 0, or alternatively, by carefully keeping track of the prefactor generated 
by the elimination of the u's and the change in the measure when going from J d[Q 
to j d[x]- The Grassmann integral for a quadratic action that is not time-translational 
invariant can in general be expressed in terms of a Pfaffian, see for example [22, 23]. 
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4.2. Consistency check 

As a consistency check we apply the formahsm of the last section to the simplest 
Majorana Hamiltonian H = —itxiX2 and compute {21x1X2) ■ In terms of the standard 
fermion annihilation operator c = ^^^"^ the Hamiltonian is H = e(c^c — 1/2) and 
2^X1X2 = 1 — 2c^c. A simple finite temperature calculation then gives (1 — 2c^c) = 
tanh(e/3/2). Let us now reproduce this result using the Majorana path integral. Taking 
derivatives of the generating function and setting the source fields to zero we obtain 

(2^X1X2) = ^) ' (53) 

i 

where in this case 

S{x, 0) = Y Xi,-u,X2,u, cos^(wj/2)-i Y Y sin(t^j)Xa,-^,Xa,a;, , (54) 

so that 

(2^X1X2) = ?7 — A^TT Y( — (55) 

M ^ 4sm (u;j72) + cos^((X'j/2)a^ 

where a = /3e/M. To perform frequency summations we use the following standard trick: 
if f{z) is a function that is regular at the points z = e~*'^^ and sufficiently well-behaved 
at infinity we have the relation 



which is often easily evaluated using residues. The contour C should enclose all of the 
singularities of f{z) but not those of {z^'^ + 1)~^. Using this, the sum in (55) is easily 
evaluated: 

f dz 1 1 
(2.X1X2) = -2«/^^(,M + i) (,_i)2_(, + i)2(,/2)2- (57) 

The residues are located at 2:1 = (1 + a/2)/(l — a/2) and Z2 = l/^i, adding up their 
contribution we arrive at 

1 1 



(2^X1X2) 



which is the correct result. 



(^f + 1) (zf + 1; 
1 1 



tanh(/3e/2), (5J 



5. Application of the Majorana path integral to the Kondo lattice 

In this section we generalize the first (favored for small J/t) variational calculation 
of section 3 to finite temperatures in terms of a path integral representation of a 
trial density matrix. Of particular interest is how to describe the formation of the 
composite object 70 = 2z/ii/i2/i3 in the path integral language. Let us first recall the 
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finite temperature variational principle [24, 25], which states that the trial free energy 
Ff, is an upper bound to the true free energy F for every properly normalized density 
matrix p, with 

Fp = TTHp + T Trplogp > F. (59) 

In our case we parametrize p through its matrix elements between coherent states in 
the rotated basis 

{^m\p\^i) OC e"^"=i^''[7lAf(k)7a,i(k)+A:,Aj(k)A„,iW]g-ELo'5'7a^ 

-^7. = E k)7:(c^„ k)7a(c^., k), (60) 

j k 

Gl\u;„k) = - (1 - ^ - /3i?^.(k)/M). 

The variational parameters are ak, which define the operators 7a(k) and /i(j(k) via (13), 
and the trial energies in the propagators E;y^{k). 

5.1. Energy expectation value 

The energy expectation values of 7's are easily calculated using the formula 

Trn.p = 1 - (^^ = " (6^) 

where a is a short-hand for (7a, k). This formula is obtained from the path integral 
using the free fermion propagator of (60), with the aid of (56) in the limit M — )■ 00. 
Using this we find 

3 

Tr i^(^)p = If, + I J]'6,(k) + X:E'^~7(k)n,.(k), (62) 



2- 

k a=l k 



and 



Ti H^^^p = - j^[Y1 sin(ak) n^,{k) - \]) [Y1 sin(ak) n^A^) - ^ 



2 

k k 

+ cyclic permutations 1 — )■ 2 — > 3 — > 1. (63) 

To arrive at these results we have also performed the average involving p's. We will now 
describe how this can be done. 

5.2. Bound state formation and the six-fermion term 

The average of terms involving p's can be formulated as integrals of the form 

J|cipi/ri)cip2/r,)rfp3/r,)j e-^^oif(;iQ) = j d[p]e-^^oH{fio). (64) 

Here -ff (po) is a short-hand for a term in the Hamiltonian where the operators pa have 
been substituted with the corresponding Grassmann number of the zeroth time slice 
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^afli cf. (41). The Grassmann numbers in (60) can be expressed in terms of the real 
space and imaginary time as 



7o%(k) = ^E-"'"""'"^^o.(r.), (65) 

loA^i) = 2^/ii,/(ri)/22,/(ri)/i3,;(ri). 

This form of the action imphes that we are gluing together the rotated Majorana 
fermions locally in space and imaginary time. To perform the integral we first make 
use the Grassmann Gaussian integral identity (i.e., a fermionic Hubbard-Stratonovich 
transformation) 

= ^ f ^[^]eE,E^^(k)G^o('^i'k)^(k)+c^.(k)7o,^^,{k)+7o%/k)^(k)]^^gg^ 



det[G'^oj J 

We then write the coupling between the auxiliary Grassmann variables ^ and the 70 as 
a sum of terms that are local in space and imaginary time 

EE'[^4(k)7o,.^,(k) +7o%(k)e.,(k)] = E[^r(r.) -ez(r.)]7o,.(r.). (67) 

j k i,l 

Before performing the integrals over ft we use another Hubbard-Stratonovich identity 
to rewrite 

^ s=± 

Here s is an auxiliary field that has an interpretation in terms of a fluctuating local 
magnetic field along the third axis. Now we use this identity at each point in space and 
imaginary time, so that we will have an S;(rj) at each point. It is also convenient to 
perform a local Z2 gauge transformation to absorb the sign into ji^ i{v.j) which will lead 
to an extra sign in the measure given by the product of all signs Y\i ■ s;(rj). We are now 
in the position to perform the integral over fii and /i2- Let us expand H^fio) in position 
and time. The most important six-fermion terms coming from (8) are 

i,d j,k,l j',k',l' 

xflifl{rj + Xd)/i2,o(rfc + Xrf)/i3,o(ri + Xd)/ii,o(r/)/i2,o(rfc')/^3,o(r/')- (69) 
We note that because of the structure of the transformation with A and B, the 
coordinates of the first three Grassmann numbers can not coincide with those of the 
last three, cf. the discussion below (18). Performing the integral over all the ftis and 
/i2's and summing over all auxiliary fields s/(rj), we see that the result is nonzero only 
when all of the coordinates coincide, i.e., when j = k = I and j' = k' = I'. For each 
point in space and time we then get a factor of 2, which can be interpreted as coming 
from the spin degeneracy. After these manipulations the coupling between ^ and ft^ in 
the exponent is simply 

E[C(r.) -e.(r.)]^3,Kr,) = ^ E'tQ (k)/i3,^^ (k) + /i^. (k)^., (k)]. (70) 

i,l j k 
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Integrating out ^ and ^* we are then back to a simple action for /is 

-^70 — ^ 5Z 5Z "^To^^^i' k), (71) 

j k 

and the 6-fermion term in question in the Hamiltonian can be replaced by 

if(6)(/io) ^^t5^5^^?,^?;/i3,o(r, + rd)/i3,o(rO. (72) 

i,d j,l 

Performing the remaining integral over J d[fLs] we generate the first term of (22), the 
second term is generated from the rotation with six i3's obtained from the original H^^'' of 
(8). The average of all other terms in the rotated Hamiltonian involving /i's will vanish. 
Physically this is due to the entirely local spin correlations of (19) that is implicit in the 
trial density matrix. 

Finally we would like to point out that the result of (71)-(72) can also be derived 
by considering the fact that the Grassmann integral picks up the coefficient of the 
term in the integrand where all Grassmann numbers are present. The expansion of 

always generates terms where the coordinates of the /t's go together, since it is an 
expansion of terms like 7o ;(ri) = 2ifL^ i{ri)fi2 i{ri)fl^ i{ri). This immediately implies that 
the coordinates of H^^\flo) has to be paired up accordingly for the result to be non- 
zero. The virtue of the formalism of this section is that it can be used also in other 
contexts where the Hamiltonian is not so simple. The presented Hubbard-Stratonovich 
decoupling schemes will be useful in e.g. Quantum Monte Carlo calculations. 

5.3. Entropy calculation 

We would now like to calculate the entropy Sp = —Tr pin p for our trial density matrix, 
which is parametrized in terms of (60). A typical term in the action is, with the short- 
hand notation a = (7a,k), 

j k,l 

Integrating out the intermediate variables {■Ca fc}fc=i,...,M-i and {C,ai}i=2,...,M we have 

This implies that the density matrix only depends on the average occupations numbers 
Ua of the states 

^^G„(a;,)e-^-^^l-n„. (75) 
j 

With the free fermion propagator of (60) the frequency sum in the limit M — )■ oo 
gives the standard parametrization Ua = {1 + e^^")'^, which also satisfy the necessary 
constraint that < fla < 1. This implies that 
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where in the last step we have used the properties of fermion coherent states, see e.g. 
[23]. Normahzing the density matrix properly it is now straightforward to calculate Sp, 
and the result for the trial entropy becomes 

3 

Sp = N\n2-Y,Y!{n% (k) ln[n^„ (k)] + [1 - n^^ (k)] ln[l - (k)] } . (77) 

a=0 k 

The first term is due to the spin system, which is maximally disordered. The second 
term is the usual entropy term for partially occupied fermion states. Extremizing the 
trial free energy we find that the trial energies should be chosen such that 

which are the standard mean field equations at finite temperature [26]. It remains to 
choose the optimal ak, which proceeds as in section 3. 



6. Conclusions and Outlook 



We have studied the Kondo lattice starting from a faithful formulation of the model in 
terms of Majorana fermions [5]. This formulation suggests a novel way of looking at the 
Kondo lattice: for strong coupling a good trial state is obtained by having "deconfined" 
independent Majorana fermions, while for weak coupling it is energetically favorable for 
the Majoranas coming from the /-spins to form a three-body bound state. The real 
situation is of course something in between, and a mix of the two pictures. In this paper 
we have investigated this scenario through a relatively simple variational calculation. 
We have also formulated a discrete imaginary time path integral formalism for Majorana 
fermions, the result (51) is both simple and beautiful. It will hopefully be useful for 
further investigations of this and other models. The investigation presented here can 
be extended in several ways. One way is to allow for non-trivial spin correlations in 
the trial state as discussed in section 3.3. Another way is to try to extend the trial 
density matrix calculation to allow for a frequency-dependent self-energy. Finally, the 
Majorana formalism could be applied to study the Kondo lattice model also away from 
half-filling, which is a straightforward extension of this work. 
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